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We construct a simple, exactly solvable model illustrating the conjecture that the evolution of 
the wave function of a system of very many interacting particles is likely to drive the one-particle 
spectrum towards a "thermal" shape. We speculate that this pseudo-thermalization, operating at the 
level of multiparton wave functions could perhaps explain the apparently thermal spectra observed 
in those high-energy processes, where the conditions for a genuine thermalization are not met. 
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Consider a multiparticle state in a finite box. If the 
state were belonging to a classical statistical ensemble, 
it would thermalizc in the presence of an interaction be- 
tween particles. In particular, the one-particle momen- 
tum distribution would tend to the Boltzmann form with 
time going on. Suppose now that the state is a pure quan- 
tum state. We wish to address the following questions: 
is it likely that the one-particle spectrum is approach- 
ing, for time — > +oo, a limit independent of the initial 
state of the system, like in a generic Markov process? 
Is it plausible that this limit distribution looks thermal, 
reflecting the features of the corresponding classical en- 
semble? Although these questions may sound academic, 
they are in fact motivated by phenomenology. However, 
since the answers are, we believe, of intrinsic theoretical 
interest, we postpone the discussion of the phenomeno- 
logical matters to the final paragraphs and focus first on 
what appears to be, hopefully, more than an amusing 
speculation. 

It is difficult to approach the problem in its full general- 
ity. Consequently, we adopt the one-particle approxima- 
tion: we follow a single test particle moving in an exter- 
nal mean field, whose action represents the multiparticle 
environement. We know, that in a classical gas this mo- 
tion generically leads to thermalization. A time arrow is 
attached to the motion of the test particle. This time 
arrow is also present when one goes over to the quantum 
description. This means that the evolution generator de- 
scribing the motion of the test particle is non-hermitian. 
Hence, it will, in general, have complex eigenvalues. Con- 
sequently, a single quantum state of the test particle, in- 
dependent of its initial state, will be selected for time 

► ~r"00. This state, to be regarded as the stationary one, 
is reached by a process, which corresponds to the decay 
of a quantum fluctuation and which is likely to resemble 
thermalization. That is the general idea of this paper. It 
will be illustrated by a simple example. 
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Consider a gas of classical particles first. How does 
the thermalization occur? A given particle of the gas is 
slowed down (accelerated) when it collides with a particle 
less (more) energetic than itself. Eventually, its momen- 
tum fluctuates near the average momentum measured by 
the gas temperature. To mock this behavior assume that 
our test particle gets random kicks, such that the transi- 
tion p — > p ± S occurs with probability 

Prob{p->p±5) ~e A±(p) (1) 

where 

A±(p) = ±[E(p)-E(p±6)] (2) 

Here E(p) is the energy of a particle. The transition 
probability satisfies the detailed balance equation 

Prob(p)Prob(p — > p ± 5) = Prob(p ± S)Prob(p ± S — » p) 

(3) 

with Prob(p) ~ exp[—E(p)/T], which is therefore the 
limit distribution of the process defined by © Equa- 
tion ([I]) violates the time reversal. It must be so, the 
process has the time arrow built into it, since it mocks 
the effect of the interaction of the test particle with its 
multiparticle environement. 

In order to simplify writing we assume first that 5 is 
a constant vector. Thus the motion is one-dimensional. 
The extension of the discussion to a spectrum of momen- 
tum transfers is trivial. 

Eq. ((T|) implies the following master equation for the 
momentum probability distribution w(p, t): 

dw(p,t)/dt = w(p + S,t)e A ^ p+5 ^ 
+ w(p - S,t)e A+(p - S) - w(p,t)[e A ~ (p) + e A+(p) ] (4) 

where the time t is measured in appropriate units, so 
that the interaction strength does not appear explicitly. 
It can be immediately seen that w(p) — e~ Eijp ^ T is a 
stationary solution of (HJ), as expected. 
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We now go over to the quantum world and replace the 
particle by a wave. We assume that the kicks are due to 
the interaction 

V = 5>-e A ~ (p+5 Ha p+5 + V + e A+ ^ala p ^ s ] (5) 
v 

where rf^ are some phase factors, scalar functions of p 
and S, while aj, (a p ) are time independent creation (an- 
nihilation) operators. The interaction operator is not 
hermitian, the time reversal invariance is broken like in 
the classical model. 

The number of particles is - by assumption - conserved 
Q • Hence the state at time t has necessarily the form 

|t>=X>(P,*)4l°> (6) 
p 

The norm of this state is time dependent, the evolution 
being non-unitary. This does not matter for our pur- 
poses, since we are interested in this paper in the shape 
of the momentum distribution only. The Schroedinger 
equation describing the evolution of the wave packet ^ 
reads 

idip{p, t)/dt = 7fe A ~ {p+5) i;{p + 5,t) 

+ r l +e A+ ( p - 5 ^(p-5,t) (7) 

In order to determine the behavior of the solution we 
have to solve first the eigenvalue equation: 

ife A ~ {p+s) ^{p + S)+ V + e A+{p ' 5) ilj{p - 6) = \4>{p) (8) 

The transition probabilities are determined by the values 
of . In the absence of a better insight into the problem 
it is reasonable to assume that rf 1 also depends primarily 
on this parameter. In order to obtain an analytically 
solvable model we postulate a linear behavior 

Arg(7 7 ± ) =e + eiA ± (pT5) (9) 
Set iP(p) = e -BM(i+^i)/2T z ( p ) t0 get 

e l6o [z(p + S) + z{p-S)]= Xz{p) (10) 

From here on the calculation becomes elementary. As- 
suming periodic boundary conditions one has 

2e^° cos (x ■ S) z(x) = \z(x) (11) 

where 

z (x) = j24py p - x (12) 

p 

Hence the sought eigenvalues and eigenvectors are 

\ x = 2e lio cos (x ■ S) (13) 

and 

lpW(p) = e -S( P )(l+2iSx)/2T CQS {p . x) (w) 

^i a) (P) = e - B (pK 1+2 ^)/ 2T sin (p • x) (15) 



When one considers a spectrum of momentum transfers, 
one gets instead of (fT3| : 

\ x = J^(<5) cos (x • <5) (16) 

a 

with some fi(5), which we assume to be a positive weight 
so that, generically, x = corresponds to the eigenvalue 
with the largest modulus. The general solution of the 
Schroedinger equation is 

V>(p,t) = 5>^i a) b) + S x 4 s Hp)]e- ltx * (17) 

X 

Let < £o < Then only the contribution corre- 
sponding to x — survives in the limit t — > +oo, be- 
cause this is the one whose norm is " explosing" the most 
rapidly (remember that the state is not normalized; the 
x = cigenstate should be regarded as the station- 
ary solution). Consequently, the asymptotic spectrum 
is w (p) - e- E ^l T . This is the generic behavior, which 
can be avoided by choosing very unnatural initial condi- 
tions only. At this point the goal we have set for our- 
selves at the beginning of this paper has been achieved. 
A rather robust quantum model has been constructed il- 
lustrating the following conjecture: the evolution of the 
wave function of a system of very many interacting par- 
ticles is likely to drive the one-particle spectrum towards 
a "thermal" shape. This behavior, which we would like 
to call pseudo-thermalization, is likely but not manda- 
tory, in the model it occurs for a (wide) range of model 
parameters only. 

Notice, that the physical contents of the classical and of 
the quantum model are different. In the classical model a 
particle keeps changing its momentum. In the quantum 
model the components of a wave keep splitting. All these 
components are a priori present at any time, only their 
amplitudes get modified. However, in both models the 
observable one-particle spectrum has the same " thermal" 
shape. Moreover, both models have the Markov prop- 
erty: the generic limiting behavior is independent of the 
initial state of the system. 

As long as one stays within the one-particle framework, 
several specific assumptions of our quantum model can be 
relaxed. Still, even in this framework, it would be inter- 
esting to check how generic is the result obtained in this 
paper, viz. does it hold for processes which are signifi- 
cantly different from |T]). The difficulty is mainly techni- 
cal, at least if one looks for an analytic solution. However, 
going beyond the one-particle picture and checking the 
conjecture in a genuine multiparticle context requires a 
significant conceptual advance. On the theory side, we 
made just a first step and much remains to be done. 

As mentioned at the beginning of this text, our re- 
search has actually a phenomenological motivation. In- 
deed, in high energy processes involving hadrons the 
inclusive spectra often look thermal, in circumstances 
where genuine thermalization seems very unlikely. For 
example, in soft hadronic collisions the inclusive trans- 
verse momentum spectrum at not too large p± has an 
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approximately boltzmannian shape, 



da 

d 2 p± 



3 -B(px)/T 



(18) 



as if it resulted from a thermalization process occuring 
during the collision. Here, E(j>±) denotes the "transverse 
energy" of a secondary [3[ and T is a "temperature" pa- 
rameter, roughly independent of the secondary hadron 
species. However, the thermalization in the final state, 
debatable but plausible in a heavy-ion collision, is nearly 
impossible in a minimum bias nucleon-nucleon collision, 
where (IT8t is the best observed. 

As is well known, in the infinite momentum frame 
P — > oo a hadron can be regarded as a 2D nonrelativis- 
tic, coherent multiparton state enclosed in a box, with 
dimensions fixed by the confinement scale. The kinetic 
energy of a parton is p]_/2xP, x denoting its longitu- 
dinal momentum fraction. It is plausible that pseudo- 
thermalization holds for the multiparton wave function. 
If so, then one can argue, invoking the parton-hadron du- 
ality, that after the hadron is broken by the interaction 
the inclusive spectrum reflects the " thermal" distribution 
holding at the parton level [f|. 

Notice, that since xP plays here the role of "mass" , one 
expects that the average transverse momentum is larger 
for larger longitudinal momentum (in cms). Such a cor- 
relation is indeed observed. On the other hand it is not 
expected when the Schwinger mechanism is operating. 

Another example: in ref. [8j a surprisingly good and 
economic description of structure functions is presented, 



using a statistical ansatz. This might be another example 
of pseudo-thermalization, this time oberved directly at 
the parton level. 

On the phenomenology side, our exercise supports the 
idea that the "thermal" spectra observed in high-energy 
processes involving hadrons do actually reflect pseudo- 
thermalization affecting the structure of the multiparton 
wave functions of the incoming hadrons. If true, this is 
to large extent independent of the details of the micro- 
scopic dynamics - a common feature of many multipar- 
ticle production phenomena - except that the validity of 
the partonpicture requires an asymptotically free theory, 
of course Q ■ 
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Note added: Preparing the original version of this 
note I have completedy overlooked results obtained in 
the context of so-called quantum chaos theory, and in 
particular the early paper by M. Srednicki [lflj, much 
more profound that the heuristic discussion presented 
above. Working with a gas of hard spheres, he has given 
strong arguments in favor of the conjecture that one- 
particle distribution in a stationary pure quantum state 
is Boltzmann, provided the corresponding classical state 
is chaotic. 
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